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ABSTRACT 

Context. Previously unremarkable, the extragalactic radio source GB 1310+487 showed a y-ray flare on 2009 November 18, reaching 
a daily flux of ~ 10 -6 photons cm -2 s -1 at energies E > 100 MeV and became one of the brightest GeV sources for about two weeks. 

Its optical spectrum shows strong forbidden-line emission while lacking broad permitted lines, which is not typical for a blazar. 
Instead, the spectrum resembles those of narrow emission-line galaxies. 

Aims. We investigate changes in the object’s radio-to-GeV spectral energy distribution (SED) during and after the prominent y-ray 
flare with the aim of determining the nature of the object and of constraining the origin of the variable high-energy emission. 

Methods. The data collected by the Fermi and AGILE satellites at y-ray energies; Swift at X-ray and ultraviolet (UV); the Kanata, 

NOT, and Keck telescopes at optical; OAGH and WISE at infrared (IR); and IRAM 30 m, OVRO 40 m, Effelsherg 100 m, RATAN-600, 
and VLBA at radio are analyzed together to trace the SED evolution on timescales of months. 

Results. The y-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift^ = 0.638. It shines through an unrelated 
foreground galaxy atz = 0.500. The AGN light is probably amplified by gravitational lensing. The AGN SED shows a two-humped 
structure typical of blazars and y-ray-loud narrow-line Seyfert 1 galaxies, with the high-energy (inverse-Compton) emission domi- 
nating by more than an order of magnitude over the low-energy (synchrotron) emission during y-ray flares. The difference between 
the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the y-ray flux and 
spectral index, with the hardest spectrum observed during the brightest y-ray state. The y-ray flares occurred before and during a slow 
rising trend in the radio, hut no direct association between y-ray and radio flares could be established. 

Conclusions. If the y-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed 
GeV spectral variability may result from varying relative contributions of these two emission components. This explanation fits the 
observed changes in the overall IR to y-ray SED. 

Key words, quasars: individual: GB 1310+487 - galaxies: jets - gamma rays: galaxies - radiation mechanisms: nonthermal - 
galaxies: active 


1. Introduction 

Blazars are active galactic nuclei (AGNs) in which relativisti- 
cally beamed emission from the jet dominates the radiative out- 
put across most of the electromagnetic spectrum. The spectral 
energy distribution (SED) of a blazar has two broad compo- 
nents: one peaking between far-IR and X-ray wavelengths and 
the other peaking at y-rays (Abdo et al. 2010b). The low-energy 
emission component is believed to be dominated by synchrotron 
radiation of relativistic electrons/positrons in the jet. Radiation 
at higher energies could be due to the inverse-Compton scat- 
tering of synchrotron photons emitted by the electrons them- 
selves (synchrotron self-Compton process, SSC; e.g., Jones et al. 
1974; Ghisellini & Maraschi 1989; Marscher & Travis 1996) 
and/or photons from external sources (external Compton pro- 
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cess, EC; e.g., Sikora et al. 1994; Dermer & Schlickeiser 2002). 
The sources of the external seed photons for the EC process 
include the accretion disk, broad-line region (BLR) clouds, 
warm dust (dusty torus), synchrotron emission from other 
faster/slower regions of the jet, and the cosmic microwave back- 
ground (CMB), with their relative contributions varying for dif- 
ferent objects. The models based on inverse-Compton scatter- 
ing by relativistic electrons are generally referred to as lep- 
tonic models (Celotti & Ghisellini 2008; Ghisellini & Tavecchio 
2009; Boettcher 2010, 2012). An alternative view regarding the 
origin of blazar high-energy emission is represented by hadronic 
models (Miicke & Protheroe 2001; Miicke et al. 2003; Sikora 
2011), where relativistic protons in the jet are the primary ac- 
celerated particles. We adopt the leptonic models as the basis for 
the following discussion. 
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Two types of radio-loud AGNs give rise to the blazar 
phenomenon: flat-spectrum radio quasars (FSRQs) and 
BL Lacertae-type objects (BL Lacs). Flat-spectrum radio 
quasars are characterized by high luminosities, prominent 
broad emission lines in their optical spectra, and the peak of 
synchrotron jet emission occurring at mid- or far-IR wave- 
lengths. Thermal emission, probably originating in the accretion 
disk surrounding the central black hole, may contribute a 
significant fraction of the optical and UV emission in some 
FSRQs (Villata et al. 2006; Jolley et al. 2009; Abdo et al. 
2010a). BL Lacertae-type objects, on the other hand, show 
mostly featureless optical spectra dominated by the nonthermal 
continuum produced by a relativistic jet. Their synchrotron 
emission peak is located between far-IR and hard-X-ray energies 
(Padovani & Giommi 1995; Fossati et al. 1998; Ghisellini et al. 
1998). In GeV y-rays, BL Lacs show a wide distribution of 
spectral slopes, while FSRQs almost exclusively exhibit soft 
y-ray spectra (Abdo et al. 2010c). It is not clear whether there is 
a physical distinction between BL Lacs and FSRQs, or if they 
represent two extremes of a continuous distribution of AGN 
properties such as black hole mass ( M .), spin, or accretion rate 
(Ghisellini et al. 2011). Recently, five radio-loud narrow-line 
Seyfert 1 galaxies (NLSyls) have been detected in y-rays by 
FermifLAT, suggesting the presence of a new class of y-ray- 
emitting AGNs (Abdo et al. 2009d; D’Ammando et al. 2012). 
The relationship between NLSyl and blazars is under debate. It 
has been suggested that radio-loud NLSyl galaxies harbor rela- 
tivistic jets (Foschini 2013; D’Ammando et al. 2013), but unlike 
blazars they are powered by less massive black holes hosted by 
spiral galaxies (Yuan et al. 2008; Komberg & Ermash 2013). 
The presence of a relativistic jet is supported by observation of 
superluminal motions in the parsec-scale radio jet of the NLSyl 
SBS 0846+513 (D’Ammando et al. 2012). The observational 
evidence that radio-loud NLSyl have M. smaller than those of 
blazars has recently been challenged by Calderone et al. (2013). 
Some nearby radio galaxies including Cen A (NGC 5128), Per A 
(NGC 1275, 3C84), and Vir A (M87, 3C274) are detected by 
FermifLAT (Abdo et al. 20 1 Oj). While part of their y-ray 
luminosity is attributed to inverse-Compton scattering of CMB 
photons on the extended (kpc-scale) radio lobes of the galaxies 
(Cheung 2007; Abdo et al. 2010i), contribution from the core 
region is also evident (Abdo et al. 2009b, 2010h). Unlike other 
radio galaxies studied by Abdo et al. (20 1 Oj), Per A exhibits 
episodes of rapid GeV variability (Donato etal. 2010; Ciprini 
2013). The core y-ray emission in radio galaxies is probably 
produced by the same mechanisms as in blazars, but with less 
extreme relativistic beaming. 

Since early satellite observations established the association 
of some discrete y-ray sources with AGNs, it became clear that 
blazars emit a considerable fraction of their total energy output 
above 100 MeV (Swanenburg et al. 1978; Hartman etal. 1999; 
Mukherjee 2002). The current generation of space-based y-ray 
telescopes that use solid-state (silicon) detectors is represented 
by instruments onboard AGILE (Tavani et al. 2009, 2008) and 
Fermi (Atwood et al. 2009), which open a window into the world 
of GeV variability and spectral behavior of y-ray-loud AGNs. In 
contrast to previous expectations (Vercellone et al. 2004), most 
of the brightest y-ray blazars detected by Fermi and AGILE were 
already known from the EGRET era (Tavani 2011). On the other 
hand, many blazars previously unknown as y-ray emitters were 
observed to reach high fluxes (> 10 -6 photons cm -2 s -1 at ener- 
gies E > 100 MeV) for only a short period of time during a flare. 
In this work, we present a detailed investigation of one such ob- 
ject. 


The radio source GB 1310+487 [also known as 
GB6B1310+4844 1 , and CGRaBS J13 12+4828, listed in 
the Fermi y-ray source catalogues as 1FGLJ13 12.4+4827 
(Abdo etal. 2010d) and 2FGLJ1312. 8+4828 (Nolan et al. 
2012), radio VLBI position 2 a J20 oo = 13 h 12 m 43 s .353644 ± 
0.22 mas, dnooo = +48 o 28'30794047 ± 0.16mas (Beasley et al. 
2002)] is a flat-spectrum radio source. It was unremarkable 
among other faint y-ray detected blazars (the E > 100 MeV 
flux during the first 11 months of the Fermi mission was 
~ 3 x 10 -8 photons cm -2 s -1 , as reported in the 1FGL catalogue; 
Abdo et al. 201 Od) until it appeared in the daily Fermi sky with 
a flux of ~ 1.0 X 10~ 6 photons cnT 2 s _1 on 2009 November 18 3 
(Sokolovsky et al. 2009). AGILE observations reported two days 
later confirmed the high-flux state of the source (Bulgarelli et al. 
2009). Follow-up observations in the near-lR (Carrasco et al. 
2009) and optical (Itoh et al. 2009) also found GB 1310+487 in 
a high state compared to historical records. The daily average 
y-ray flux remained at ~ 1.0 X 10 -6 photons cm -2 s _1 for more 
than a week (Hays & Escande 2009). 

This paper presents multiwavelength observations of 
GB 1310+487 before, during, and after its active y-ray state, and 
suggests possible interpretations of the observed SED evolution. 
In Sect. 2 we describe the observing techniques and data anal- 
ysis. Sect. 3 presents an overview of the observational results. 
In Sect. 4 we discuss their implications, and we summarize our 
findings in Sect. 5. Throughout this paper, we adopt the fol- 
lowing convention: the spectral index a is defined through the 
energy flux density as a function of frequency F v oc y +Q ', the 
photon index r p h is defined through the number of incoming 
photons as a function of energy dN(E)/dE oc E r +, and the 
two indices are related by T p h = 1 - a. We use a ACDM cos- 
mology, with the following values for the cosmological param- 
eters: Hq = 71kms~ 1 Mpc~*, F2 m = 0.27, and fi A = 0.73 (see 
Komatsu et al. 2009; Hogg 1999), which corresponds to a lumi- 
nosity distance of Di = 3800 Mpc, an angular-size distance of 
Da = 1400 Mpc, and a linear scale of 6.9 pc mas -1 at the source 
redshiftz = 0.638 (see Sect. 3.4). 

2. Multiwavelength observations 

2. 1. Gamma-ray observations with Fermi/LAT 

Fermi Gamma-ray Space Telescope (FGST; hereafter Fermi) 
is an orbiting observatory launched on 2008 June 11 by a 
Delta II rocket from the Cape Canaveral Air Force Station in 
Florida, USA. The main instrument aboard Fermi is the Large 
Area Telescope (LAT; Atwood et al. 2009; Abdo et al. 2009a; 
Ackermann et al. 2012), a pair-conversion telescope designed to 
cover the energy band from 20 MeV to greater than 300 GeV. 
The Fermi/ LAT is providing a unique combination of high sen- 
sitivity and a wide field of view of about 60°. Fermi is operated 
in an all-sky survey mode most of the time, which makes it ideal 
for monitoring AGN variability. 

The dataset reported here was collected during the first 
33 months of Fermi science observations from 2008 August 4 
to 2011 June 13 in the energy range 100 MeV - 100 GeV. The 
33 month time interval is divided into subintervals according to 
the level of its y-ray activity as observed by FennifLAT (see Ta- 
ble 1). 

1 The correct B1950 source name, if its declination is expressed as 
three digits, is 1310+487, while 1310+484 is an unrelated nearby radio 
source. 

2 see http://astrogeo.org/vlbi/solutions/rfc_2012b 

3 UT dates are used through the text. 
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Table 1. Changes in the y-ray spectrum between time intervals considered in the analysis. 

Period UT interval 100 MeV-100 GeV flux TS N - 

33 months 2008-08-04-2011-06-13 (1.03 + 0.04) x 10” 7 2.18 + 0.02 44T5 3566“ 

pre-flare 2008-08-04-2009-11-16 (0.34 ± 0.05) x 10“ 7 2.41 +0.09 205 476 

Flare 1 2009-11-16-2009-12-21 (6.94 + 0.32) X 10“ 7 1.97 + 0.03 3333 947 

Interflare 2009-12-21-2010-04-26 (1.37 + 0.11) x 10“ 7 2.15 + 0.06 917 592 

Flare 2 2010-04-26-2010-07-26 (2.83 + 0.16) X 10“ 7 2.14 + 0.04 1839 907 

post-flare 2010-07-26-2011-06-13 (0.44 ± 0.06) x 10~ 7 2.34 + 0.09 236 422 

Column designation: Col. 1, y-ray activity state; Col. 2, time interval used for spectral analysis; Col. 3, average flux in units of 
photons cm -2 s -1 ; Col. 4, photon index: dN(E)/dE oc E~ r p h ; Col. 5, Test Statistic (TS) defined in Sect. 2.1; and Col. 6, number of 
photons attributed to the source (model dependent). 


FermifLAJ data comprise a database containing arrival 
times, directions, and energies of individual silicon-tracker 
events supplemented by information about the spacecraft posi- 
tion and attitude needed to calculate the effective exposure for 
a celestial region and time interval of interest. The maximum- 
likelihood method is used to analyze these data by constructing 
an optimal model of the sky region as a combination of point- 
like and diffuse sources having a spectrum associated with each 
one of them (Mattox et al. 1996; Abdo et al. 20 1 Od). The sig- 
nificance of source detection is quantified by the Test Statis- 
tic (TS) value, determined by taking twice the logarithm of the 
likelihood ratio between the models including the target source 
(Z-i ) and one including only the background sources (Lq): TS = 
2(lnZj -InZo). The ratios Lq and L\ are maximized with respect 
to the free parameters in the models. The Monte-Carlo simula- 
tion performed by Mattox et al. (1996) for EGRET confirmed 
theoretical predictions (Wilks 1938) that for a GeV telescope, in 
most cases, the TS distribution is close to^ 2 . 

The unbinned likelihood analysis was performed with 
the Fermi Science Tools package 4 version v9r21p0. The 
DIFFUSE class events in the energy range 100 MeV - 100 GeV 
were extracted from a region of interest defined as a circle 
of 15° radius centered at the radio position of GB 1310+487. 
A cut on zenith angle > 100° was applied to reduce 

contamination from Earth-limb y-rays, produced by cosmic 
rays interacting with the upper atmosphere (Shawetal. 2003; 
Abdo et al. 2009c). Observatory rocking angles of greater than 
52° were also excluded. A set of instrument response func- 
tions (IRFs) P6_V11_DIFFUSE was used in the analysis. The 
sky model contained point sources from the 2FGL catalogue 
(Nolan et al. 2012) within 20° from the target, as well as 
Galactic gll_iem_v02_P6_Vll_DIFFUSE . fit and isotropic 
isotropic_iem_v02_P6_Vll_DIFFUSE.txt diffuse compo- 
nents 5 . All point-source spectra were modeled with a power law; 
the photon index was fixed to the catalogue value for all sources 
except the target. The diffuse-background parameters were not 
fixed. The estimated systematic uncertainty of flux measure- 
ments with LAT using P6 V11 IRFs is 10% at 100 MeV, 5% 
at 500 MeV, and 20% at 10 GeV and above 6 . 

The lightcurve of the target source was constructed by ap- 
plying the above analysis technique to a number of indepen- 


dent time bins. The time bin width was chosen to be seven days. 
Sources with less than one photon detected in the individual bin 
or with TS < 25 were excluded from the sky model for that 
bin. The lightcurves were computed by integrating the power- 
law model in the energy range 100 MeV - 100 GeV. 

For lightcurves with time bins of fixed widths, the choice 
of bin width is a compromise between temporal resolution and 
signal-to-noise ratio for the individual bins. For Fermi / LAT 
an alternative method has recently been developed (Lott et al. 
2012), in which the time bin widths are flexible and chosen to 
produce bins with constant flux uncertainty. Flux estimates are 
still produced with the standard LAT analysis tools. In this case 
we used Fermi Science Tools v9r27pl and P7CLEAN_V6 
event selection and IRFs, for which the current version of the 
adaptive binning method has been optimized (we have checked 
that using the P7S0URCE_V6 class yields very similar fluxes). At 
times of high source flux, the time bins are narrower than dur- 
ing lower flux levels, therefore allowing us to study more rapid 
variability during these periods. 

The lower energy limit of the integral fluxes computed for 
the adaptively binned lightcurve is chosen to minimize the bin 
widths needed to reach the desired relative flux uncertainty for 
most bins. The derivation of this energy limit, called the op- 
timum energy, is presented by Lott et al. (2012). Because the 
source is variable and the optimum energy value depends on 
the flux, we compute the optimum energy with the average flux 
over the first two years of LAT operation reported in the 2FGL 
catalogue (Nolan et al. 2012). The optimum energy is found to 
be E 0 - 283 MeV for this source. We produced two sets of 
adaptively binned lightcurves in the 283 MeV - 200 GeV energy 
range, one with 25% flux uncertainties and another with 15% 
uncertainties. For each of these uncertainty levels we created a 
second version of the lightcurve by performing the adaptive bin- 
ning in the reverse-time direction. 

2.2. Gamma-ray observations with AGILE/GRID 

The AGILE y - ray satellite (Tavani et al. 2009, 2008) was 
launched on 2007 April 23 by a PSLV rocket from the Satish 
Dhawan Space Centre at Sriharikota, India. AGILE is a mis- 
sion of the Italian Space Agency (ASI) devoted to high-energy 


— ; astrophysics, and is currently the only space mission capa- 

For documentation of the Science Tools, see k le 0 f observing cosmic sources simultaneously in the en- 
K ://f ^' 9SfC '^ S ,a ■f v/s f/d ata/ a nalys is /docl me ntat i o a / bands 18 _ 60ke y and 30MeV _ 30GeV thanks t0 its 

3 The models are available from the Fermi Science Support Center x . .- . . , 

,,r ■ r //i*./ /-i i/D i two scientific instruments: the hard X-ray Imager (Super- 
http ://fermi .gsfc.nasa.gov/ssc/data/access/lat/BackgroundiJoflels.frtmT . r,nnn\ j r- n t • V, . 

6 For newer VI y 6 IRFs used in the adaptive lightcurve analysis AGlre ; Peroci et a 1 2007) and the Gamma-Ray Imaging Detec- 

described below, the systematic uncertainties are lower at high ener- tor (GRID, Rappoldi & AGILE Collaboration 2009). During the 

gies: 10% at 100 MeV, 5% at 560 MeV, 10% at 10 GeV and above first two years of the mission, AGILE was mainly operated by 
(Ackermann et al. 2012); however, this difference is not critical for the performing 2-4 week-long pointed observations, but following 
present analysis. the reaction wheel malfunction in October 2009 it was operated 
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in a spinning observing mode, surveying a large fraction of the 
sky each day. 

The AGILE/ GRID instrument detected enhanced y- ray 
emission from GB 1310+487 from 2009 November 20 17:00 
(JD 2455156.2) to 2009 November 22 17:00 (JD 2455158.2) 
(see Bulgarelli et al. 2009, for preliminary results). Level 1 
AGILE/ GRID data were reanalyzed using the AGILE Standard 
Analysis Pipeline (see Pittori et al. 2009; Vercellone et al. 2010, 
for a description of the AGILE data reduction). We used y- 
ray events from the ASDCSTDe archive, filtered by means of 
the FM3.119 pipeline. Counts, exposure, and Galactic back- 
ground y-ray maps were created with a bin size of 0!5 X 0°5 , for 
E > 100 Me V. Since AGILE was in its spinning observing mode, 
all maps were generated including all events collected up to 50° 
off-axis. We rejected all y- ray events whose reconstructed di- 
rections form angles with the satellite -Earth vector smaller than 
90°, reducing the y-ray Earth limb contamination by excluding 
regions within ~ 20° from the Earth limb. We used the latest ver- 
sion (BUILD-2®) of the Calibration files (IS02 3), which will be 
publicly available at the ASI Science Data Centre (ASDC) site 7 , 
and the y-ray diffuse emission model (Giuliani et al. 2004). We 
subsequently ran th e AGILE Multi-Source Maximum Likelihood 
Analysis (ALIKE) task using a radius of analysis of 10° in order 
to obtain the position and the flux of the source. A power-law 
spectrum with a photon index T = 2.1 was assumed in the anal- 
ysis. 

2.3. X-ray observations with Swift/XRT 

The X-ray Telescope (XRT; Burrows et al. 2005) onboard the 
Swift satellite (Gehrels et al. 2004) provides simultaneous imag- 
ing and spectroscopic capability over the 0.2-10 keV energy 
range. The source GB 1310+487 was observed by Swift at seven 
epochs during the two y-ray activity periods and in June 2011 
during the low post-flare state. A summary of the Swift ob- 
servations is presented in Table 2. Swift/XRT was operated in 
photon-counting (pc) mode during all observations. The low 
count rate of the source allows us to neglect the pile-up effect 
which is of concern for the XRT in pc mode if the count rate 8 is 
> 0.6 count s _1 . 

The xrtpipeline task from the HEASoft v6 . 14 package 
was used for the data processing with the standard filtering 
criteria. To increase the number of counts for spectral analy- 
sis, the resulting event files were combined with xselect to 
produce average X-ray spectra for the periods of Flare 1 and 
Flare 2 defined in Table 1 . The spectrum for the Flare 1 pe- 
riod was binned to contain at least 25 counts per bin to utilize 
the x 2 minimization technique. The combined spectra were an- 
alyzed with XSPEC vl2 .8.1. The simple absorbed power-law 
model with the H I column density fixed to the Galactic value 
Ah i = 0.9 1 7 X 10 20 cm -2 (obtained from radio 21 cm measure- 
ments by Kalberla et al. 2005) provided a statistically acceptable 
fit (reduced^ 2 = 1.2 for 12 degrees of freedom) to the 0.3- 
lOkeV spectrum. Leaving Ah i free to vary results in the values 
Ah i = 0.4-0. 5 x 10 22 cm -2 and r p h x-ray ~ 1.3. However, this 
model does not improve the fits. The low photon counts prevent 
a more detailed study. 

Individual observations obtained during the periods of 
Flare 1 and Flare 2 were also analyzed using the same fixed-An i 
model, but no evidence of spectral variability within the periods 
was found; however, the low photon counts could easily hide 

7 http://agile.asdc.asi.it 

s http : //www. swift . ac.uk/pileup . shtml 


mild spectral changes. The Swift /X RT observation obtained dur- 
ing the Flare 2 and post-flare intervals (Table 1) resulted in a low 
number of detected photons. The Cash (1979) statistic is applied 
to fit this dataset with the absorbed power-law model. The Cash 
statistic is based on a likelihood ratio test and is widely used for 
parameter estimation in photon-counting experiments. The net 
count rate in the 0.3-10 keV energy range changed by a factor of 
1.6 between the two observations conducted during Flare 1 and 
by a factor of 3.5 over the whole 33 -month period. The X-ray 
spectral analysis results are presented in Table 2. 

2.4. Uitraviolet-optical observations 

The Swift Ultraviolet-Optical Telescope (UVOT; Roming et al. 
2005) has a diameter of 0.3 m and is equipped with 
a microchannel-plate intensified CCD detector operated in 
photon-counting mode. Swift/ UVOT observed GB 1310+487 si- 
multaneously with Swift/XRT. Various filters were used at dif- 
ferent epochs ranging from the U to M2 bands (as detailed 
in Table 2), with the best coverage achieved in the U band. 
Since the object is very faint, multiple subexposures taken during 
each observation were stacked together with the tool uvotimsum 
from the HEASoft package. A custom-made script based on 
uvotsource was employed for aperture photometry (using the 
standard 5" aperture diameter) and count rate to magnitude con- 
version taking into account the coincidence loss (pile-up) correc- 
tion (Poole et al. 2008; Breeveld et al. 2010). The Galactic red- 
dening in the direction of this source is E(B - V) - 0.013 mag 
(Schlegel et al. 1998). Using the extinction law of Cardelli et al. 
(1989) and coefficients presented by Roming et al. (2009), the 
following extinction values were obtained for the individual 
bands: A v = 0.041, A B = 0.053, A v = 0.065, and A M i = 0.122 
mag. Magnitude-to-flux-density conversion was performed us- 
ing the calibration of Poole et al. (2008). 

A star-like object is visible in Nordic Optical Telescope 
(NOT) images about 3" southwest of the AGN. This object 
would be blended with the AGN in UVOT images which lack 
adequate angular resolution. Contribution of this object to the 
total flux measured by UVOT is the likely reason for the dis- 
crepancy between UVOT and NOT U- band measurements dur- 
ing the low state of GB 1310+487. The nearby galaxy (Sect. 3) 
also contributes to the measured UVOT flux. 

The Nordic Optical Telescope, a 2.5 m instrument located 
on La Palma, Canary Islands, conducted photometric observa- 
tions of GB 1310+487 with its ALFOSC camera on 2010 July 7 
and 2011 May 29 during the second flare and the post-flare low 
state, respectively. The VaST 9 software (Sokolovsky & Lebedev 
2005) was applied for the basic reduction (bias removal, flat- 
fielding) and aperture photometry of the NOT images. A fixed 
aperture 175 in diameter was used for the measurements. 
The source 3UC 277-1 16569, which served as the compari- 
son star for the Kanata observations (see below), was satu- 
rated on NOT i band images and so could not be used. Instead, 
SDSSJ131240. 83+482842. 9 (a J2 ooo = 13 h 12 m 40 s .84, d J2 ooo = 
+48°28'4279, Abazajian et al. 2009; see Fig. 1) was used as 
the comparison star. Its Johnson-Cousins magnitudes were com- 
puted from the SDSS photometry using conversion formulas of 
Jordi et al. (2006): U = 19.133 ± 0.081,5 = 19.303 ± 0.013, 
V = 18.822 + 0.012,5 = 18.589 + 0.01 1, and 1 = 18.177 + 0.018 
mag. 

Kanata, a 1.5 m telescope at the Higashi-Hiroshima Ob- 
servatory, observed GB 1310+487 in the R and I bands with 

9 http://scan.sai.msu.ru/vast 
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Table 2. Swift observations of GB 1310+487. 


ObsID 

Date (UT) 

JD 

Exposure 

UVOT 

XRT 

0.3-10 keV unabs. flux 

Tph X-ray 



2455... 

(ks) 

(mag) 

(cts/s) 

(10 -13 erg cm -2 s _1 ) 

Flare 1 period 

001 2009-11-27 

002 2009-11-30 

162.77 

165.97 

8.7 

8.1 

M2 =21.1(2) 

U = 21.2(4), B> 21.3, V > 20.9 

0.022(2) 

0.013(1) 

16.7 ±3.3 

0.88 ±0.18 

Flare 2 period 
003 2010-06-25 

372.66 

4.5 

U > 21.3, B > 21.5, V > 20.5 

0.010(1) 



004 

2010-07-03 

380.65 

4.5 

U = 20.9(3), B > 21.5, V > 20.3 

0.009(1) 

6.6+ 1.5 

1.15 ±0.23 

005 

2010-07-07 

384.63 

4.9 

U = 21.4(3) 

0.009(1) 

006 

2010-07-11 

388.68 

4.6 

(7 = 21.3(3) 

0.010(1) 



post-flare period 
007 2011-06-03 

717.12 

9.5 

U = 21.1(2) 

0.006(1) 

5.5 + 1.8 

0.93 ±0.34 


Column designation: Col. 1, observation number in the Swift archive omitting the leading 00031547; Cols. 2, 3, date of observation 
given by the Gregorian and Julian Date; Col. 4, exposure time in kiloseconds; Col. 5, SwiftfUWOT photometry (here and later in 
the text the error in parentheses corresponds to the last decimal place of the value before the parentheses); Col. 6, SwiftpLKT net 
count rate in counts/s and its uncertainty; Col. 7, 0.3-10 keV unabsorbed flux derived from fitting Swift/XKT data with the power- 
law model (datasets 1-2 and 3-6 are combined to increase the photon statistics); and Col. 8, X-ray spectral index (r p h x-ray) as 
defined in Sect. 1 . 


o 


5 " 


Fig. 1 . Nordic Optical Telescope 7?-band image of the GB 1310+487 
region obtained on 2011 May 29. The exposure time was 300 s. 
North is up and east is to the left. The AGN and comparison star 
SDSSJ13 1240.83+482842.9 used in the NOT data analysis are marked 
with the letters “Q” and “C,” respectively. 


the HOWPol instrument (Kawabata et al. 2008) in the nonpo- 
larimetric mode for eight nights during the first and second 
y - ray flares. Relative point-spread function (PSF) photometry 

was conducted using 3UC 277-1 16569 (orj 2 ooo = 13 h 12 m 54 s .09, 

<5j2ooo = +48°27'5872, J2000; R = 16.109, 1 = 15.657 
mag; Zacharias et al. 2010) as the comparison star. The adopted 
Galactic extinction values were Ar = 0.035 and A / = 0.025 mag 
(Schlegel et al. 1998). The calibration by Bessell et al. (1998) 
was employed for the magnitude-to-flux conversion. 

The source GB 1310+487 was assigned a redshift of 0.501 
based on a 2007 March 21 1200 s Hobby-Eberly Telescope 
Low Resolution Spectrograph (HET/LRS) observation 10 , which 
showed strong [O II] 43727 at 5592 A and weak evidence of host 
absorption features (Healey et al. 2008; Shaw et al. 2012). This 
spectrum had insufficient signal-to-noise ratio (S/N) to exclude 
weak broad lines, or to cleanly measure the optical continuum, 
leaving the nature of the source uncertain. 

Thus, we reobserved the source with the Keck 10 m tele- 
scopes. Long-slit spectra were obtained with the Keck II DEep 


10 Falco et al. (1998) previously reported:: = 0.313, but it was indicated 
as a “marginal measurement.” 


Imaging Multi-Object Spectrograph (DEIMOS) (Faber et al. 
2003) on 2013 April 07 with ~ 1" seeing. Two 600 s inte- 
grations were obtained using the 600 lines per mm (7500 A 
blaze 11 ) grating, providing coverage in the 4450-9635 A range 
with a ~ 100 A gap between the two CCDs. With the 170 slit, 
the spectra have an effective resolution of ~ 3.0 A. Conditions 
were good, but not completely photometric; the flux scale might 
be uncertain by roughly a factor of 2. Moreover, we obtained 
2 x 180 s g- and /Gbaud images of the object with the two cam- 
eras on the Keck I Low Resolution Imaging Spectrometer (LRIS) 
(Oke et al. 1995) on May 10 under ~ 1" seeing, as shown in Fig- 
ure 2. The DEIMOS slit on April 07 was placed on the bright 
core of the source, at the parallactic angle (Filippenko 1982) of 
PA = 143° (measured from N to E), and the extended wings of 
the host were also included. The companion was ~ 3" off the 
slit. 

A second Keckll/DEIMOS spectrum was obtained on 
June 10 with a different slit position and ~ 0.7" seeing. It has 
a higher signal-to-noise ratio than the first DEIMOS spectrum; 
however, it was affected by a cosmic-ray hit that prevented accu- 
rate measurement of H-y in the z = 0.638 system, and conditions 
were not photometric when the standard star was being observed. 
The two spectra are normalized to epoch 1 (April 07) using the 
[O II] 43727 line atz = 0.500. The continuum cannot be used to 
cross-calibrate the two spectra because of the significantly vari- 
able AGN flux contribution; the second-epoch continuum level 
appears to have dropped relative to the emission lines by ~ 1/3. 
The two spectra were averaged for further analysis. 

2.5. Infrared photometry 

Observations in the near-IR were carried out with the 2.1m tele- 
scope of the Guillermo Haro Observatory, INAOE, Mexico. The 
telescope is equipped with the CANICA camera together with J, 
H, and K s filters. We carried out differential photometry between 
the object of interest and other objects in the 5' x 5' field. The 
observations showed an increase of about one magnitude dur- 
ing the Flare 2 period with respect to the Flare 1 and post-flare 
periods (results are summarized in Table 3). Magnitudes are re- 


11 The “blaze wavelength” is the wavelength for which the grating is 
the most efficient. 


Article number, page 5 of 18 



A&A proofs: manuscript no. 1310 





Fig. 2. Keckl/LRIS images of GB 1310+487 obtained on 2013 May 10. 
Top panel: the g- and R-band images (8" X 8" field of view; north is 
up and east is to the left). Bottom panel: g- and R -hand images after 
subtraction of a scaled point source from the offset core. The boxes 
show the location of the DEIMOS slit in the 2013 April 07 observations. 
The circles indicate the location of the radio-loud AGN. These residual 
images reveal a relatively undisturbed foreground galaxy. 


ferred to the 2MASS 12 survey published photometry. The source 
GB 1310+487 itself was not detected in the 2MASS survey. The 
survey detection limit is J - 15.8, H — 15.1, and K s = 14.3 mag 
(Skrutskie et al. 2006). These values may be considered upper 
limits on the object brightness at the 2MASS observation epoch 
ofJD 2451248.8408 (1999 March 11). 

The source GB 1310+487 is listed in the Wide-field In- 
frared Survey Explorer (WISE; Wright et al. 2010) catalogue 
(Cutrietal. 2012) with the following magnitudes in the four 
WISE bands: 3.4/rm W1 = 12.302 ± 0.024, 4.6/an W2 = 
11.254 + 0.021, 12 pm W3 = 8.596 + 0.021, and 22 pm. W4 
= 6.368 + 0.044. The IR colors (W1-W2 = 1.048 + 0.032, W2- 
W3 = 2.658 + 0.030 mag) are at the blue edge of the area in 
the color-color diagram occupied by blazars and Seyfert galax- 
ies (see Fig. 12 in Wright et al. 2010, Fig. 1 in D’Abrasco et al. 
2012, and Fig. 1 in Massaro et al. 2011), indicating that the AGN 
and not the host galaxy’s stars or warm dust is responsible for 
most of the IR flux in these bands. WISE observations of this 
area were conducted on 2010 June 3-8 during Flare 2 (Table 1). 

2.6. Radio observations 

As part of an ongoing blazar monitoring program, the Owens 
Valley Radio Observatory (OVRO) 40 m radio telescope has 
observed GB 1310+487 at 15 GHz regularly since the end of 
2007 (Richards et al. 2011). This monitoring program studies 
over 1500 known and likely y-ray-loud blazars, including all 
CGRaBS (Healey et al. 2008) sources north of declination -20°. 
The objects in this program are observed in total intensity twice 
per week. The minimum measurement uncertainty is 4mJy 
while the typical uncertainty is 3% of the measured flux. Obser- 

12 http : //www. ipac . caltech . edu/2mass/ 


Table 3. Ground-based photometry of GB 1310+487. 


Date 

JD (UTC) 
2455... 

Filter 

mag 

Instrument 

Flare 1 period 
2009-11-28 

164.33668 

R 

20.58(2) 

Kanata 

2009-11-29 

165.33620 

R 

20.86(9) 

Kanata 

2009-12-05 

171.31613 

R 

20.61(8) 

Kanata 

2009-11-28 

164.34442 

I 

19.41(2) 

Kanata 

2009-11-29 

165.35134 

I 

19.66(6) 

Kanata 

2009-12-05 

171.31613 

I 

>18.91 

Kanata 

2009-12-13 

179.33910 

I 

19.52(1) 

Kanata 

2009-11-22 

158.04390 

H 

15.87(6) 

OAGH 

interflare period 
2010-03-17 272.95922 

H 

15.97(5) 

OAGH 

Flare 2 period 
2010-07-07 

385.44392 

U 

21.91(6) 

NOT 

2010-07-07 

385.47325 

B 

22.6(1) 

NOT 

2010-07-07 

385.43551 

V 

21.62(5) 

NOT 

2010-06-03 

351.03843 

R 

>20.78 

Kanata 

2010-06-04 

352.05203 

R 

>20.52 

Kanata 

2010-06-05 

353.02920 

R 

>20.64 

Kanata 

2010-07-07 

385.47773 

R 

20.85(2) 

NOT 

2010-07-17 

395.07192 

R 

>19.90 

Kanata 

2010-07-19 

396.99981 

R 

>20.48 

Kanata 

2010-06-03 

351.04981 

I 

19.98(7) 

Kanata 

2010-06-05 

353.04057 

I 

19.80(2) 

Kanata 

2010-07-07 

385.46087 

I 

19.79(2) 

NOT 

2010-07-17 

395.08708 

I 

>19.27 

Kanata 

2010-06-15 

362.70970 

J 

16.4(1) 

OAGH 

2010-06-16 

363.75472 

J 

16.6(1) 

OAGH 

2010-06-18 

365.69752 

J 

16.3(1) 

OAGH 

2010-06-19 

366.71342 

J 

15.7(1) 

OAGH 

2010-05-17 

333.79362 

H 

14.59(1) 

OAGH 

2010-05-20 

336.83470 

H 

15.11(5) 

OAGH 

2010-06-15 

362.70186 

H 

14.61(7) 

OAGH 

2010-06-16 

363.74752 

H 

14.68(3) 

OAGH 

2010-06-19 

366.69745 

H 

14.84(5) 

OAGH 

2010-06-15 

362.72145 

K s 

13.8(1) 

OAGH 

2010-06-16 

363.76233 

K s 

13.78(8) 

OAGH 

2010-06-19 

366.72800 

K s 

13.7(1) 

OAGH 

post-flare period 
2011-05-29 711.46544 

V 

21.56(4) 

NOT 

2011-05-29 

711.46135 

R 

20.79(2) 

NOT 

2011-05-29 

711.46922 

I 

19.87(3) 

NOT 

2011-07-31 

773.69450 

H 

15.83(7) 

OAGH 


Column designation: Cols. 1, 2, the Gregorian and Julian Date 
of observation, respectively; Col. 3, filter; Col. 4, magnitude and 
its uncertainty; and Col. 5, telescope name. 


vations are performed with a dual-beam (each 2.5' full width at 
half-maximum intensity, FWHM) Dicke-switched system using 
cold sky in the off-source beam as the reference. Additionally, 
the source is switched between beams to reduce atmospheric 
variations. The absolute flux-density scale is calibrated using ob- 
servations of 3C286, adopting the flux density (3.44 Jy) from 
Baars et al. (1977). This results in a ~ 5% absolute flux-density- 
scale uncertainty, which is not reflected in the plotted errors. 

Multifrequency radio observations of GB 1310+487 were 
performed with the 100 m Effelsberg telescope operated by the 
MPIfR 13 . Observations were conducted on 2009 December 1 


13 Max-Planck-Institut fur Radioastronomie 
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following the reported y- ray flare, on 2010 June 28 during the 
second y-ray active phase, and on 2011 June 5. Secondary fo- 
cus heterodyne receivers operating at 2.64, 4.85, 8.35, 10.45, 
and 14.60GHz were used. The observations were conducted 
with cross-scans (i.e., the telescope’s response was measured 
while slewing over the source position in azimuth and eleva- 
tion). The measurements were corrected for (a) pointing off- 
sets, (b) atmospheric opacity, and (c) elevation-dependent gain 
(see Fuhrmann et al. 2008; Angelakis et al. 2009). The multifre- 
quency observations were completed within 1 hr for each observ- 
ing session. The absolute flux-density calibration was done by 
observing standard calibrators such as, 3C48, 3C 161, 3C286, 
3C 295, and NGC 7027 (Baars et al. 1977; Ott et al. 1994). 

The 1RAM 30 m Pico Veleta telescope observations at 86.24 
and 142.33 GHz took place on 2009 December 7. The observa- 
tions and data-reduction strategy were similar to those with Ef- 
felsberg; a detailed description is given by Nestoras et al. (2014). 
Both the Effelsberg 100 m and IRAM 30 m telescope observa- 
tions were conducted in the framework of the F-GAMMA pro- 
gram (Fuhrmann et al. 2007; Angelakis et al. 2008, 2010, 2012; 
Fuhrmann et al. 2014). 

For comparison with the latest Effelsberg, IRAM, and 
OVRO results, we use data from the RATAN-600 576 m 
ring radio telescope of the Special Astrophysical Observatory 
(Russian Academy of Sciences); RATAN-600 observations of 
GB 1310+487 were performed in transit mode at the southern 
sector with the flat reflector quasi-simultaneously at 3.9, 7.7, 

11.1, and 21.7 GHz in June 2003 within the framework of the 
spectral survey conducted by Kovalev et al. (1999b, 2002). The 
flux-density scale is set using calibrators listed by Baars et al. 
(1977); Ott et al. (1994). Details on the RATAN-600 obser- 
vations and data processing are discussed by Kovalev et al. 
(1999a). 

The National Radio Astronomy Observatory’s Very Long 
Baseline Array (VLBA, Napier 1994, 1995) is a system of ten 
25 m radio telescopes dedicated to very long baseline interfer- 
ometry (VLBI) observations for astrophysics, astrometry, and 
geodesy. After publication of the report on the November 2009 
y- ray flare (Sokolovsky et al. 2009), GB 1310+487 was added 
to the MOJAVE 14 program (Lister et al. 2009a). Three epochs of 
VLBA observations at 15 GHz were obtained between 2009 and 
2010 . 

3. Results 

3. 1. y-ray analysis 

The y-ray counterpart of GB 1310+487 was localized, inte- 
grating 33 months of Fermi/LXT monitoring data, to aj 2 ooo = 
198° 187, fcooo = 48/472, with a 68% uncertainty of 0/014 = 
50". This is a factor of five larger than the spacecraft align- 
ment accuracy of 10" = 0/003 (Nolan et al. 2012). The y-ray 
position is only 0/005 = 18" away from the radio position 
of GB 1310+487. Within the Fermi/LAT error circle no other 
radio sources are seen with the VLA FIRST 1.4 GHz survey 
(White et al. 1997), which provides the best combination of sen- 
sitivity and angular resolution for that region of the radio sky 
to date. Therefore, the positional association of the y-ray source 
with the radio source GB 1310+487 is firmly established. The X- 
ray brightening observed during the first and, to a lesser extent, 
the second y-ray flares together with the near-IR brightening dur- 
ing the second y-ray flare support the identification of the y-ray 

14 Monitoring Of Jets in Active galactic nuclei with VLBA Experi- 
ments, http : / /www . physics . purdue . edu/astro/MO J AVE/ 



Fig. 3. Fermi/LAT y-ray SEDs compared to the power-law models 
(shown as lines in this logarithmic plot) derived from unbinned like- 
lihood analysis. We stress that the power-law models are not fits to the 
binned energy flux values; rather, these are two independent ways to 
represent Fermi/LAT photon data. The time intervals defined in Table 1 
are color labeled. 


source with the lower-frequency counterpart. See Tables 1, 2, 3, 
and the discussion below for details. 

The y-ray spectra of GB 1310+487 at various activity states 
listed in Table 1 are presented in Figure 3. The plotted spectral 
bins satisfy the following requirements: TS > 50 and/or model- 
predicted number of source photons N > 8. The flux value at 
each bin was computed by fitting the model with source position 
and power-law photon index fixed to the values estimated over 
the entire period in the 0.1-1 00 GeV energy range. The Galactic 
component parameters were fixed, as were all other nontarget 
source components of the model. 

To test if the power law (PL) is an adequate approxima- 
tion of the observed y-ray spectrum, the PL fit (represented 
by a straight line if plotted on a logarithmic scale) was com- 
pared to the fit with a log-parabola (LP) function defined as 
dN/dE = /Vo(F/Fo) _( “ +/31n(£/£o) ' , where N is the number of pho- 
tons with energy E, No is the normalization coefficient, Eq is 
a reference energy, a is the spectral slope at energy E o, and // 
is the curvature parameter around the peak. For the combined 
33 month Fermi/ LAT dataset, the fit with the assumption of a 
PL spectrum for the target source provides its detection with 
the Test Statistic TSpl = 4415, while the LP spectrum leads 
to TSlp = 4423. These values may be compared by defining, 
following Nolan et al. (2012) and in analogy with the source 
detection TS described in Sect. 2.1, the curvature Test Statistic 
TS cu rve = 2(lnZ, L p-lnZ PL ) = TSlp-TSpl. The obtained value of 
TS curve = 8 corresponds to a 2.8cr difference, which is lower than 
the TS cur ve > 16 (4 cr) threshold applied by Nolan et al. (2012). 
We conclude that while there is a hint of spectral curvature, it 
cannot be considered significant. 

The broken power-law (BPL) model was also tested, but it 
did not provide a statistically significant improvement over the 
PL or the LP models (TSbpl = 4423, for the best-fit break energy 
E b = 3 GeV, the photon indexes T p hi = 2.30 ± 0.04, r p hi = 
0.05 + 0.02 above and below the break, respectively). Therefore, 
we adopt the simpler PL model for the following analysis. 

The 33 -month y-ray lightcurve of the source obtained with 
the seven-day binning is presented in Figure 4. Two major flar- 
ing periods are clearly visible. The first, brighter flare peaked 
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around 2009 November 27 (JD 2455 163) with the weekly av- 
eraged flux of (1.4 ± 0.1) X 10~ 6 photons cm~ 2 s _1 . The peak 
flux averaged over the two-day interval centered on that date 
is (1.9 ± 0.2) x 10 -6 photons cm -2 s -1 . The source continued to 
be observed at a daily flux of ~ 0.5 X 10~ 6 photons cm~ 2 s _1 
for another two weeks. The second flare peaked around 2010 
June 17 (JD 2455365) at the seven-day integrated flux of 
(0.54 ± 0.07) X 10 -6 photons cm -2 s -1 . The daily flux of ~ 
0.5 X 10~ 6 photonscnV 2 s _1 was observed for about three weeks 
around this date. The two flares demonstrate remarkably con- 
trasting flux evolution: the first is characterized by a fast rise 
and slower decay, while the second flare shows a gradual flux 
rise followed by a sharp decay. Following Burbidge et al. (1974), 
Valtaoja et al. (1999), and Gorshkov et al. (2008), we define the 
flux-variability timescale as t v ar = At/AlnA, where A In S is 
the difference in logarithm of the photon flux at two epochs 
separated by the time interval At. The observed flux-variability 
timescale during the onset of Flare 1, as estimated from the 
seven-day binned lightcurve (Fig. 4), is t mT =» 3 days. The 
timescale of flux decay after Flare 2 is f var ~ 5 days. 

The Fermi/ LAT lightcurve constructed with the alternative 
analysis method, the adaptive binning (with 25% flux uncer- 
tainty at each bin), is presented in Figure 5. It confirms all the 
features visible in the constant bin-width lightcurve, but also al- 
lows us to investigate fast variability during high-flux states in 
greater detail. The first flare episode, Flare 1, consists of four 
prominent subflares, each with a time width of a day or less. 
The subflares show no obvious asymmetry and the variability 
timescale f var for the three point rise of the second and third 
subflares (at JD 2455157.5 and 2455161.5) was estimated as 
0.36 ± 0.20 days and 0.45 + 0.23 days, respectively. A second 
adaptively binned lightcurve was produced in the reverse-time 
direction, which gives a similar, but not identical, time binning. 
The result of the timescale estimates for this second version of 
the lightcurve was found to be consistent with the first analy- 
sis. A similar analysis for the lightcurves with 15% uncertainties 
give timescale estimates of about 1 day for the most rapid vari- 
ability. We conclude that the adaptively binned lightcurves show 
evidence of a variability timescale of half a day with a conserva- 
tive upper limit of 1 day. For the second and fainter flare epoch 
the timescales seen in the adaptive binning are consistent with 
the estimate from the fixed-binned lightcurve described above. 

Table 1 presents spectral analysis results for the different y- 
ray activity states of the source: “pre -flare” and “post-flare” peri- 
ods represent the low-activity level, “Flare 1” and “Flare 2” rep- 
resent the high-activity state, while during the “interflare” inter- 
val the source showed an intermediate y- ray flux level. Figure 3 
presents the observed Fermi [LAX spectrum during these states. 
Significant evolution of the y-ray photon index, T p h, is detected 
between the different flux states (Table 1). Figure 6 presents T p h 
as a function of ( E > lOOMeV) flux. The harder- when-brighter 
trend is clearly visible. 

Integrating the AGILE observations from 2009 Novem- 
ber 18 12:00 (JD 2455154.0) to 2009 November 22 12:00 UT 
(JD 2455158.0), we obtain a y-ray flux 7 7 £>iooMeV = (1.65 ± 
0.48) x 1 0 -6 photons cm -2 s -1 , at a significance of VTS = 6.1. 
This result is in good agreement with the flux value derived 
from the preliminary analysis by Bulgarelli et al. (2009). Prior 
to the Fermi launch, AGILE observed GB 1310+487 (in pointing 
mode) during two other periods, but did not detect the source. 
During the first period (from 2007 October 24 12:00 UT to 2007 
November 1 12:00 UT, JD 2454398.0-2454406.0), the 2cr upper 
limit was F e> ioo MeV < 0.28 X 10 -6 photons cm -2 s _1 , while in 
the second period (from 2008 April 30 12:00 UT to 2008 May 10 
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Fig. 4. Weekly binned FermifLPG lightcurve. Blue filled circles are val- 
ues with TS > 25, gray filled circles are values with 5 < TS < 25, and 
gray arrows indicate 2cr upper limits for time bins with no significant 
detections (TS < 5). A four-day integrated AGILE data point is added 
as an open box for comparison. 
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JD-2450000.0 

Fig. 5. Fermi/LAA lightcurve constructed with the adaptive binning 
method (Lott et al. 2012). The magnified plot of Flare 1 is shown in the 
insert. The energy range for this lightcurve is chosen to minimize the 
uncertainties in time and flux, while the lightcurve was Fig. 4 is given 
in the commonly used E > 100 MeV energy range. 


12:00 UT, JD 2454587.0-2454597.0) we obtained a 2cr upper 
limit of Fe>\oo MeV < 0.31 X 10~ 6 photons cm~ 2 s -1 . 

3.2. X-ray to infrared spectrum 

Results of the X-ray spectral analysis are presented in Table 2. 
The obtained values of the X-ray photon index r p h x-ray are 
among the hardest reported for blazars (Giommi et al. 2002; 
Donato et al. 2005; Sikora et al. 2009). Radio-loud NLSyl have 
r p h x-ray similar to the ones found in blazars (Paliya et al. 2013; 
Abdo et al. 2009d). However, it cannot be excluded that the X- 
ray spectrum with an intrinsic value of r p h X -ray is artificially 
hardened by additional absorbing material along the line of sight 
(see the discussion of NOT imaging results below). Future high- 
quality X-ray observations are necessary for investigating this 
possibility. 
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Fig. 6. The y-ray photon index, T P h, as a function of flux for the time 
periods defined in Table 1 . 


The UV and blue parts of the optical spectrum are flat, in 
contrast to the steep spectrum seen in the near-IR ( iJHK bands). 
Also, the observed variability amplitude is decreasing toward 
bluer wavelengths, with the [7-band brightness being essentially 
constant. 

The H- band flux showed an increase of about one magnitude 
during the Flare 2 period with respect to the Flare 1 and post- 
flare periods, in contrast to the behavior seen in other bands. The 
results of ground-based photometric measurements are summa- 
rized in Table 3. 


3.3. Imaging with NOT 

The Nordic Optical Telescope images (Fig. 1) show a fuzzy 
extended object, probably a galaxy, with a point source offset 
076 from its center. Considering that there are many galaxies 
of comparable brightness visible in the field, this picture may 
be interpreted as the AGN (corresponding to the point source) 
shining through an unrelated foreground galaxy. This may be 
the source of confusion in the AGN’s redshift determination 
(Sokolovsky et al. 2009; Flealeyetal. 2008; Falcoetal. 1998), 
and it also explains the steepness of the optical-lR SED (the ob- 
served SED was corrected for Milky Way absorption, but ab- 
sorption in the intervening galaxy may also be significant). On 
the other hand, it is not uncommon for AGN host galaxies to 
have disturbed morphologies, making it appear that the AGN is 
off-center. 

The galaxy contributes a large fraction of the total optical 
flux, when the point source is in the low state. If the host galaxy 
of GB 1310+487 is similar to the giant ellipticals studied by 
Sbarufatti et al. (2005), (Mr) = -22.9 ± 0.5 mag, its magnitude 
at z = 0.500 should be R ~ 20 (or 0.6 mag fainter at z = 0.638; 
Sect. 3.4). A typical NLSyl from the Veron-Cetty & Veron 
(2010) catalogue having (My) = -21 .4 mag would appear 1 mag 
fainter than a giant elliptical in the R band assuming V — R — 
0.5 mag (Xanthopoulos 1996). Therefore, the observed galaxy 
could be the host of GB 1310+487. The visible offset between 
the point source and the center of extended emission could re- 
sult from the disturbed morphology of the host, as noted above. 



\(A) 


Fig. 7. Keck II/DE1MOS spectrum of GB 1310+487. The two narrow- 
line systems are indicated. The region 6900-7000 A is lost to a gap 
between the two CCDs. The slit position is shown in Figure 2. 


Table 4. GB 1310+487 emission-line strengths. 


Species T Res t 


Flux; 


0.500 


51 W 0.500 


Flux; 


MIL 


EW* 

^ _ 0.638 


[Oil] 3726+3729 4.24 + 0.07 32 ± 6 0.97 + 0.06 4.2 + 0.5 

Fly 4340 0.65 + 0.053.1 +0.7 cosmic-ray hit 

H p 4861 1.69 + 0.06 6.6 + 0.8 0.24 + 0.06 4.1 ± 0.7 

[O III] 4959 not detected 1.00 + 0.06 3.6 + 0.8 

[OIII] 5007 1.11 ±0.05 6.5 + 0.6 2.61 +0.05 12.0 ±2.2 


“ Fluxes are in units of 10 _l/ ergcm _2 s , with lcr statistical 
errors. The overall flux scale is uncertain by up to a factor of 2, 
but the relative fluxes are much more accurate. b The equivalent 
width (EW) in angstroms. 


3.4. Keck imaging and spectroscopy 

Standard reductions, extractions, and calibrations of the 
DEIMOS data produced the spectrum shown in Figure 7; it is the 
average of two observations conducted on April 07 and June 10, 
2013. The strongest line, [O II] A3727, confirms the FIET red- 
shift identification at z = 0.500, and we also see [O III] and 
narrow Balmer emission for this system. However, there are ad- 
ditional lines, mostly in the red half. These represent a second 
system with narrow forbidden and Balmer emission, this time at 
z = 0.638. The line strengths are given in Table 4. The [O II] dou- 
blets are barely resolved, but the oxygen and Balmer line widths 
are consistent with the instrumental resolution. The [O III] emis- 
sion at z = 0.638 appears resolved with a deconvolved width 
of ~ 200kms _1 . Unfortunately, the red limit of the spectrum 
does not include the Ha/[N II] lines for either system. For the 
z = 0.500 system, we cover [O I] T6300, which is weak or ab- 
sent. For the z = 0.638 system, we cover Mg II T2800, and can 
place a 3 <x rest equivalent width limit of ~ 1.0 A on any broad 
emission; the H fi line is marginally detected for the z = 0.638 
system at 4<x level. The ratio of H/J to [O III] (z = 0.638) is 
small, even if the H/7 flux is treated as an upper limit. Together 
with the resolved [O III] this indicates nuclear excitation. 
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Thus, we clearly have two superimposed systems and wish 
to identify which system hosts the radio-loud core (and, by infer- 
ence, the y-ray source). The Keck I LRIS images confirm the ba- 
sic structure seen in the NOT images; the source is extended with 
a brighter core displaced ~ 076 to the west. Figure 2 shows 8" 
regions around the AGN. The DE1MOS slit position on April 07 
is marked on the g frames (left). At the bottom we show the im- 
ages after removal of a point-source PSF ( g - 23.89, R = 22.45 
mag) from the offset core. The residuals show a relatively regu- 
lar galaxy having FWFIM = 178, with g = 21.95 and R = 20.59 
mag. The coordinate system was referenced through the SDSS 
image of the field, with an estimated uncertainty relative to the 
radio frame of 072 ; the circles show the position of the VLB1 
source (Sect. 1) and have radii twice this uncertainty. Flence, the 
radio source is coincident with the point-like peak of the com- 
bined source. We also find that the z = 0.500 emission lines are 
offset 0721 + 0706 SE along the slit from the z = 0.638 system, 
toward the continuum tail representing the extended galaxy. The 
deprojected offset is ~ 0735 E of the AGN core. We thus con- 
clude that the true AGN redshift is z = 0.638, and we are view- 
ing it through an approximately face-on galaxy showing strong 
narrow-line emission. 

Our extracted spectrum is weighted toward the AGN core, 
although it also contains appreciable light from the foreground 
galaxy. Both spectra are dominated by narrow forbidden lines, 
yet there is appreciable continuum associated with both compo- 
nents as well. The foreground galaxy is probably not an AGN, 
but we cannot be certain; without the [N HJ/Hor line ratio, we 
are unable to fully distinguish “LINER” (Low Ionization Nu- 
clear Emission-line Region) emission from an H II region (e.g., 
Ho et al. 1997). However, the strong [O II] 43727 and lack of ob- 
vious [O I] 46300 argue against a power-law ionizing spectrum, 
suggesting that the z = 0.500 emission represents star formation 
in the foreground galaxy lacking AGN activity. 

The positional accuracy of the available observations of mul- 
tiwavelength variability (Sect. 2) is not sufficient to distinguish 
between the foreground and background objects discussed here 
and in Sect. 3.3 as the source of high-energy emission. The pro- 
posed interpretation that the background AGN is the high-energy 
source rests on the consideration that the observed fast y-ray 
variability (Sect. 3.1) is typical of radio-loud AGNs (which the 
background source is), while there are no firm indications of 
AGN activity in the foreground galaxy. 

3.5. Results of radio observations 

The radio spectrum of GB 1310+487 is generally flat, with a 
wide peak located between 22 GHz and 86 GHz (Table 5). The 
variability amplitude at 2.64 GHz is slightly lower compared to 
higher frequencies. The 15 GHz lightcurve of GB 1310+487 ob- 
tained with the OVRO 40 m telescope and complemented by 
measurements with the Effelsberg 100 m and the VLBA is pre- 
sented in Figure 8. It shows a period of high activity with two 
separate peaks that started in mid-2010 and is still ongoing. 

The 15 GHz VLBA images (Fig. 9) show two emission re- 
gions separated by ~ 0.4 mas. To quantify their parameters we 
fit the observed visibilities with a model consisting of two circu- 
lar Gaussian components using the Difmap software (Shepherd 
1997). The modeling results are presented in Table 6. The uncer- 
tainties in parameters of the model components were estimated 
following Lee et al. (2008), and the resolution limit achieved for 
each component was computed following Lobanov (2005) and 
Kovalev et al. (2005). 


Table 5. Multifrequency radio observations of GB 1310+487. 


V 

(GHz) 

F v 

(Jy) 

(T 

(Jy) 

V 

(GHz) 

Fv 

(Jy) 

cr 

(Jy) 

RATAN-600 

2003-06 

Effelsberg 100 m 

2010-06-28 

21.74 

0.401 

0.104 

2.64 

0.135 

0.002 

11.11 

0.222 

0.012 

4.85 

0.107 

0.001 

7.69 

0.157 

0.019 

8.35 

0.111 

0.002 

3.95 

0.221 

0.057 

10.45 

0.112 

0.005 

Effelsberg 1 00 m 

2009-12-01 

Effelsberg 100 m 

2011-06-05 

2.64 

0.161 

0.001 

2.64 

0.201 

0.006 

4.85 

0.133 

0.001 

4.85 

0.189 

0.003 

8.35 

0.130 

0.002 

8.35 

0.207 

0.005 

10.45 

0.130 

0.003 

10.45 

0.213 

0.010 

14.60 

0.121 

0.006 




IRAM 30 m 2009-12-07 




86.24 

0.282 

0.070 




142.33 

0.206 

0.065 





Column designation: Col. 1, the central observing frequency; 
and Cols. 2, 3, the observed flux density and its uncertainty, re- 
spectively. 



JD - 2450000.0 

Fig. 8. Radio lightcurve at 15 GHz obtained with the OVRO 40 m 
telescope (points) supplemented with two 14.6 GHz measurements ob- 
tained with the Effelsberg 100 m telescope (square). VLBA measure- 
ments of the core (component CO, Table 6) are indicated as diamonds. 
The two arrows mark the peaks of the y-ray flares observed by Fermi. 

The SW component increased its brightness during the three 
epochs. Comparison with the lightcurve in Figure 8 shows that 
this component is responsible for most of the flux observed with 
single-dish instruments. The fainter component located to the 
NE is gradually fading. If the SW component is the 15 GHz core, 
the position of the second component aligns nicely with the ori- 
entation of the kiloparsec-scale jet observed with the VLA at 
1.4 GHz by Machalski & Condon (1983). No significant proper 
motion could be detected between the three 15 GHz MOJAVE 
epochs. The 3<x upper limit which can be placed on proper mo- 
tion is q < 0.3masyr~', corresponding to jS ap P <11 (yS app is in 
units of the speed of light) at the source redshift, which is within 
the range of apparent jet speeds occupied by y-ray-bright blazars 
(Lister et al. 2009b; Savolainen et al. 2010). The projected linear 
size of the double structure resolved with the VLBA is ~ 2.7 pc 
= 8 X 10 18 cm. The overall 15 GHz VLBI polarization of the 
source measured by MOJAVE is 2. 8-6.5% which is indicative of 
beamed blazar emission. The weakly beamed, high viewing an- 
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Fig. 9. VLBA radio image of GB 1310+487 obtained on 2010-12-24 
at 15 GHz during the course of the MOJAVE program. The image map 
peak is 0.206 Jy beam -1 and the first contour is 0. 15 mjy beam -1 . Ad- 
jacent contour levels are separated by a factor of 2. Naturally weighted 
beam size is indicated at the lower-left corner of the image. Green cir- 
cles indicate positions and best-fit sizes of the model components pre- 
sented in Table 6. 

Table 6. Parsec-scale components observed at 15 GHz. 



Comp. 

Distance 

FWHM 

Flux density 

T b 


(mas) 

(mas) 

(Jy) 

(K) 


2009-12-26 = 

JD 2455192 


CO 


< 0.39 

0.090 + 0.021 

> 3 x 10 9 

Cl 

0.32 ± 0.11 

< 0.42 

0.014 + 0.004 

>4x 10 8 


2010-08-06 = 

JD 2455415 


CO 


< 0.36 

0.144 + 0.033 

> 6 x 10 9 

Cl 

0.36 + 0.10 

< 0.40 

0.008 + 0.002 

> 3 x 10 8 


2010-12-23 = 

JD 2455554 


CO 


< 0.40 

0.211 + 0.049 

> 7 x 10 9 

Cl 

0.43 ± 0.12 

< 0.52 

0.004 + 0.002 

> 8 x 10 7 


Column designation: Col. 1, component name, where CO is the 
presumed core and Cl is the decaying jet component; Col. 2, 
projected distance from the core (CO); Col. 3, FWHM of the 
Gaussian component; Col. 4, component flux density; and Col. 5, 
observed brightness temperature. 


gle sources in MOJAVE tend to be unpolarized (Lister & Homan 
2005). 

3.6. SED during the two flares 

The SED of GB 1310+487 is presented in Figure 10. It has the 
classical two-humped shape with the high-energy hump domi- 
nating over the synchrotron hump during the first (brighter) flare 
by a Compton dominance factor of q > 10. For the Flare 2 period 
the value of Compton dominance may be measured accurately 
thanks to simultaneous observations of Fermi/LXY and WISE: 
q = 12. The uncertainty of this measurement is limited by the 
accuracy of the absolute calibration of the two instruments and 


should be less than 10%. Fast variability within the Flare 2 pe- 
riod may also contribute to the uncertainty. The value of q is 
probably larger for Flare 1 than for Flare 2, judging from the 
lower near-IR flux observed during Flare 1. 


4. Discussion 


4.1. y-ray luminosity variability, and spectrum 

The monochromatic y-ray energy flux averaged over the du- 
ration of the first flare is vF v ~ 10 10 erg cm 2 s *. At the 
redshift of the source this corresponds to an isotropic lumi- 
nosity of ~ 10 47 ergs -1 . Considering the expected bolomet- 
ric correction of a factor of a few, the y-ray luminosity of 
GB 1310+487 is comparable to that typically observed in flar- 
ing y-ray blazars (e.g., Tanaka et al. 2011; Abdo et al. 20 10a, e) 
and NLSyl (D’Ammando et al. 2013). It is about two orders of 
magnitude lower than the most extreme GeV flares of 3C 454.3 
in November 2010 (Abdo et al. 2011b) and PKS 1622-297 in 
June 1995 (Mattox et al. 1997). The outstanding y-ray flare of 
3C 120 in November 1968 had a comparable isotropic lumi- 
nosity of ~ 10 47 ergs -1 (Volobuev et al. 1972). The exceptional 
GeV photon flux of 3C 120 was due to the relative proximity 
of the source (z = 0.033; Michel & Huchra 1988) compared to 
the brightest y-ray blazars mentioned above. The observed large 
y-ray luminosity of GB 1310+487 is an indirect indication of a 
high Doppler boosting factor of the source (Taylor et al. 2007; 
Pushkarev et al. 2009). 

The difference in lightcurve shape, overall duration, and 
shortest observed variability timescale between the two flares of 
the source may indicate that they occurred in different jet regions 
or were powered by different emission mechanisms as discussed 
below. In both cases, this implies differences in the emitting- 
plasma parameters for the two flares, such as the electron en- 
ergy distribution, magnetic field strength, bulk Lorentz factor, or 
external photon field strength. Variability timescales of 3 days 
and shorter are common in GeV blazars (e.g., Mattox et al. 
1997; Abdo et al. 2010f, 2011b; Sbarrato et al. 2011). The light- 
travel-time argument limits the y-ray emitting region size r < 
cSt vaT /( 1 + z) « a few x 10 15 cm, where c is the speed of light 
in vacuum, z is the source redshift, and conservatively assuming 
the Doppler factor 6 = [ F( I -J3 cos 0)] -1 < a few (we have no ev- 
idence of extreme Doppler boosting from VLB1 and y-ray data; 
Sect. 4.2), where T is the Lorentz factor, /J is the bulk velocity of 
the emitting blob in units of the speed of light, and 6 is the angle 
between the blob velocity and the line of sight. 

It is important to check that the observed harder-when- 
brighter trend in the y-ray spectrum is not related to the expected 
correlation between the flux and the index in the power-law 
model. If the number density of photons arriving from the source 
is dN/dE = Nq(E / E o)-^ (where No, and E o are constants), the 
integrated photon flux between energies E nrdx and E mm is 

f _ NqEq £™x Ll>+1 _ £min Lrph + 1 
_ r p h + 1 Eq Eq 


Assuming that T > 1 and E max is large, 
p, ~ NoEo Emin 

-r ph + 1 Eo 


The derivative 

II 


d F 
dr ph 


= F 


1 

— F p h + 1 


-In 



Eo 


, 0 
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Fig. 10. Quasi-simultaneous radio to y-ray SED of GB 1310+487 during the two flaring episodes and the post-flare period covered by our multi- 
wavelength observations. The time intervals corresponding to these events are defined in Table 1. 


if \n{E m i n /Eo) , l/(-Tph + 1). The range of parameters derived 
from our analysis is T = 1.97-2.41, E mm = 100 Me V, Eq = 
283 Me V, and -1.04 = In (E min /E 0 ) < l/(-r ph + 1) = -1.03 to 
-0.70, so d/'/dTph > 0. The expected correlation between d/-’ 
and T p h due to their mathematical dependence is positive, which 
is opposite to what is actually observed. We conclude that the 
observed harder-when-brighter trend is real and not related to 
the intrinsic correlation of the model parameters. 

Previously a harder-when-brighter trend (Fig. 6) has 
been seen at GeV energies only in a handful of blazars: 
3C273, PKS 1502+106, AO 0235+164, and 4C +21.35 by 
Fermi (Abdo et al. 2010k,f,g; Tanaka et al. 2011); 3C 454.3 
(Ackermann et al. 2010; Vercellone et al. 2010; Abdo et al. 
2011b; Stern & Poutanen 2011) and PKS 1510-089 by 
Fermi and AGILE (Abdo et al. 2010a; D’Ammando et al. 
2011); and 3C279 (Hartman et al. 2001) and PKS 0528 + 134 
(Mukherjee et al. 1996) by EGRET. The same harder-when- 
brighter behavior was suggested by the combined analysis of 
relative spectral index change as a function of relative flux 
change in a few of the brightest FSRQs, low-, and intermediate- 
peaked BL Lacs using the first six months of Fermi data by 
Abdo et al. (2010g). The current detection presents one of the 
clearest examples of this spectral behavior. 

The spectral evolution during the Flare l-interflare-Flare2 
periods (Fig. 3 and Table 1) may be qualitatively understood as 
the gradual decrease in energy of the y-ray emission peak. Dur- 
ing Flare 1, the spectrum is hard, implying that the spectral peak 
is located above or around 5 GeV. Consequently, during the in- 
terflare period the spectrum is softer, with a hint of curvature; 
the emission peak may be located around 1-2 GeV. Later, during 
Flare 2, the spectrum is also soft with a peak possibly located at 
even lower energies. This interpretation is inspired by the visual 


inspection of Figure 3. The peak-frequency evolution is difficult 
to quantify owing to the insufficient number of collected pho- 
tons, which results in the simple PL fit (with no curvature) be- 
ing a statistically acceptable model for the LAT data. However, 
the overall SED (Fig. 10) suggests that the high-energy emis- 
sion peak should be located somewhere around the LAT band. 
We can use the LAT spectral index (Tph) vs. Compton peak fre- 
quency (v£ akHz ) correlation log 10 v^ akHz = -4.0T ph + 31.6 re- 
ported by Abdo et al. (2010b) to estimate that vj£ ak ^ changed 

from 10 22 to 10 24 Hz between the pre-flare and Flare 1 periods. 
The change in vj^ ak Hz may result not from a continuous shift 
of a single y-ray emission peak, but from a change in relative 
strengths of two emission components peaking at different fre- 
quencies, as discussed in Sect. 4.9. 


4.2. Jet Doppler factor 


The Doppler factor of the relativistic jet in GB 1310+487 maybe 
constrained using two independent lines of argument: one based 
on the requirement that the emitting region should be transparent 
to its own y radiation (since we observe it), the other based on 
the absence of apparent proper motion seen by the VLBA. 

The minimum Doppler factor needed to avoid y - y atten- 
uation for y-rays interacting with lower energy photons present 
inside the emitting region may be calculated using Eq. (39) of 
Finke et al. (2008), 


6yy ^ 


2 fl_1 (l +zy- la cr T D 
m,.( A far 


L eif yn 


where it is assumed that the synchrotron flux is well repre- 
sented by a power law of index a (J^ yn oc e°), crj is the scat- 
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tering Thomson cross-section, Di is the luminosity distance 
to the source, m e is the electron mass, and E\ — E/(m e c 2 ) is 
the dimensionless energy of a y-ray photon with energy E for 
which the optical depth of the emitting region r yy = 1 (see also 
Dondi & Ghisellini 1995). The maximum energy of observed y- 
ray photons that can be attributed to the source is ~ lOGeV 
(Fig. 3), so e, = 10 GeV/(5.11 x 10“ 4 )GeV = 2 x 10 4 . This 
means = 5.11 x 10 -5 , and the corresponding frequency for 
this is 6.3 X 10 15 Flz. From the observed SED (Fig. 10), we esti- 
mate / sy 1 n w 10 -14 erg s _1 cm -2 and a ~ -2. Taking t var = 3 days 

f l 

(Sect. 3.1) we obtain 6 yy > 1.5. 

Assuming the angle between the jet axis and the line of sight, 
6, is 0 max , the one that maximizes the apparent speed, /? app , for a 
given intrinsic velocity, /?, we may estimat^ the corresponding 

Doppler factor ^vlba < 11 (if £vlba = T - y6 app + 1). We note 
that if 9 is smaller than 0 max , the actual 6 will be larger than the 
above estimate (see, e.g., Cohen et al. 2007; Kellermann et al. 
2007; Marscher 2009 for a discussion of relativistic kinematics 
in application to VLBI). 

Recent RadioAstron (Kardashev et al. 2013) Space-VLBI 
observations of high brightness temperatures in AGNs suggest 
that the actual jet flow speed is often higher than the jet pattern 
speed (Sokolovsky 2013). These results question the applicabil- 
ity of 6 estimates based on VLBI kinematics. The available lower 
limits on the core brightness temperature, 7),, in GB 1310+487 
(Table 6) are consistent with negligible Doppler boosting within 
the standard assumption of the equipartition inverse-Compton 
limited 7* « a few x 10 n K (Readhead 1994). 

4.3. Black hole mass 

If we equate the linear size estimated from the shortest ob- 
served variability timescale to the Schwarzschild radius, the 
corresponding black hole mass would be M, « 10 10 M o . 
However, TeV observations of ultra-fast (timescale of minutes) 
variability in blazars PKS 2155-304 (Aharonian et al. 2007; 
Abramowski etal. 2010) and Mrk501 (Albert et al. 2007) lead 
to M, estimates inconsistent with those obtained by other 
methods (Begelman et al. 2008), unless an extremely large 
Doppler factor 6 ~ 100 is assumed for the y- ray emitting re- 
gion (Ghisellini & Tavecchio 2008; Sbarrato et al. 2011). Short 
timescale variability may arise from the interaction of small (size 
r < r s ) objects such as stars (Barkov et al. 2012) or BLR clouds 
(Araudo et al. 2010) with a broad relativistic jet. This should 
caution us against putting much trust in the above M. estimate. 

4.4. UV, optical, and IR emission 

The UV-to-lR behavior of the source may be understood if the 
near-IR light is dominated by the synchrotron radiation of the 
relativistic jet, while in the optical-UV the contribution of line 
emission and/or thermal emission from the accretion disk starts 
to dominate over the synchrotron radiation. The line and thermal 
emission are not relativistically beamed and, therefore, more sta- 
ble compared to the beamed synchrotron jet emission, decreas- 
ing the variability amplitude in the parts of the SED where their 
contribution to the total light is comparable to that from the jet. 
Thermal emission features are observed in SEDs of many FSRQ- 
type blazars (e.g., Villata et al. 2006; Hagen-Thorn et al. 2009; 
Abdo et al. 2010a; D’Ammando et al. 2011). Starlight from the 
host galaxy also contributes to the total optical flux in some y- 
ray loud AGN (Nilsson et al. 2007; Abdo et al. 201 la, c). This 


contribution is significant mostly for BL Lac -type blazars and 
non-blazar AGN. Finally, as discussed above, the nearby galaxy 
and the star-like object, both probably unrelated to the source un- 
der investigation, may contribute to its total optical flux if an ob- 
servation lacks angular resolution to separate contributions from 
these objects. 

4.5. Radio properties 

The radio loudness parameter, R ra di 0 , defined as the ratio of 
5 GHz flux density, L 5 ghz, to the 5-band optical flux density, L B , 
is R radio = 7.5GHz /Lb ~ 10 4 . This is an order of magnitude larger 
than typical R ra di 0 values found in quasars (Kellermann et al. 
1989) and radio-loud NLSyl galaxies (Doi et al. 2006), but it is 
comparable to the largest observed values (Singal et al. 2013) 15 . 
The extremely low optical luminosity compared to the radio lu- 
minosity may either be an intrinsic property of this source, or it 
may result from absorption in the intervening galaxy (Sect. 3.3, 
3.4). 

The radio spectrum of GB 1310+487 (Table 5) is typical for 
a blazar. Relatively rapid (timescale of months) and coherent 
changes across the cm band suggest that most of the observed 
radio emission comes from a compact region no more than a 
few parsecs in size. Comparison of the 15 GHz lightcurve pre- 
sented in Figure 8 with the 15 GHz VLBA results (Table 6) in- 
dicates that the component CO (presumably the core) is the one 
responsible for most of the observed single-dish flux density of 
the source. Specifically, CO is the site of the major radio flare 
peaking around JD 2455500 (October-November2010). 

The presence of correlation between cm-band radio and y- 
ray emission is firmly established for large samples of blazars 
(Ackermann et al. 2011; Arshakian et al. 2012; Linford etal. 
2012; Kovalev 2009). The typical y-ray/radio time delay 
ranges from 1 month to 8 months in the observer’s frame, 
with y-rays leading radio emission (Pushkarev et al. 2010; 
Leon-Tavares et al. 2011). However, for individual sources it 
is often difficult to establish a statistically significant corre- 
lation because of the limited time span of simultaneous y- 
ray-radio data compared to a typical duration of radio flares 
(Max-Moerbeck et al. 2012). This could also limit our knowl- 
edge of the maximum possible radio/y-ray time delay. 

In the case of GB 1310+487, no clear connection is visible 
between its radio and y- ray activity, based both on the available 
single-dish (Fig. 8) and VLBI monitoring data (Table 6). 

4.6. Object classification 

Shaw et al. (2012) classified the optical spectrum of 
GB 1310+487 as a LINER, which is inconsistent with the 
y- ray and radio loudness (Sect. 4.5; Giuricin et al. 1988). The 
absence of broad lines precludes classification as a quasar. 
Prominent forbidden emission lines are not typical of BL Lac- 
type objects. Therefore, while being similar to blazars in its 
high-energy, radio, and IR properties, GB 1310+487 cannot 
be classified as a classical blazar on the basis of its optical 
spectrum. 

As discussed in Sect. 3.4, the point-source emission at z = 
0.638 observed by Keck is most likely related to the AGN. Pogge 
(2000) defines NLSyl as having permitted lines only slightly 
stronger than forbidden lines, [O III]/H/? < 3, and FWHM(H/T) 

15 Singal et al. (2013) use the rest frame luminosity at 5 GHz and 
2500 A respectively to define R. This definition of R should be con- 
sistent within a factor of a few with the one we use. 
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< 2000 km s -1 . The anomalously strong [O III] T/14959, 5007 
emission formally disqualifies this source, and tends to support 
a Seyfert 2 (or a narrow-line radio galaxy, considering the ob- 
ject’s radio loudness) classification, which would be difficult to 
understand if the radio- and y-ray jets align with the Earth’s line 
of sight. Our S/N is too low to allow unambiguous detection of 
Fe II emission. Broad H/3, if present, is weaker than a third of the 
narrow component, and there is no evidence of Mg II 2800 A. 
Thus, no broad-line component is observed. We also find that 
Ca H&K are weak, if present, and the 4000 A break is smaller 
than 0.1. These aspects suggest appreciable nonthermal lumi- 
nosity for the core AGN emission. Thus, we tentatively advance 
the view that synchrotron emission from the AGN dominates the 
variable point-source core, but that a surrounding narrow line 
region dominates the line flux. 

Giirkan et al. (2014) studied the WISE infrared colors of 
radio-loud AGN; GB 1310+487 falls in a region of the color di- 
agram occupied mostly by quasars and broad-line radio galax- 
ies, although some narrow-line radio galaxies are also present. 
The source GB 1310+487 is well away from the locus of low- 
excitation radio galaxies (LERGs) and also has a 22 pm luminos- 
ity (~ 4 X 10 45 ergs -1 ) typical of high-excitation radio galaxies 
(HERGs; see Fig. 8 in Giirkan et al. 2014). However, if one uses 
the criteria of Jackson & Rawlings (1997) GB 13 10+487 would 
qualify as a LERG based on its optical spectrum. We note that 
WISE photometry of the AGN might be contaminated by the 
foreground galaxy. 

Being a narrow-line radio-loud AGN, the object is not a 
member of common types of y-ray flaring extragalactic sources 
(blazars and NLSyls). One possibility is that the object is anal- 
ogous to nearby radio galaxies like Per A with additional ampli- 
fication due to gravitational lensing that makes y-ray emission 
from its core detectable at high redshift. The similarity to Per A 
is supported by its lack of superluminal motion (Lister et al. 
2013), low 8 inferred from SED modeling (Abdo et al. 201 Oh), 
and absence of changes in VLBI and single-dish radio properties 
that can be attributed to GeV events (Nagai et al. 2012). 

Another possibility is that the object may be a bona fide 
blazar with its optical non-thermal emission swamped by the 
host elliptical as proposed by Giommi et al. (2013) as possible 
counterparts of unassociated Fermi sources. In this case, how- 
ever, one would not understand the observed variable optical 
point source. Higher S/N spectroscopy with increased wave- 
length coverage would be helpful in characterizing the z = 0.638 
y- ray/radio AGN. Higher resolution spatial imaging is needed to 
probe the nature of the foreground (z = 0.500) galaxy. 


4.7. Gravitational lensing 


Considering that the AGN is located behind the visible disk of 
another galaxy (Sect. 3.4, 4.6), amplification of the AGN light by 
gravitational lensing is a real possibility. In the simplified case 
of a point lens, the AGN light is amplified by a factor of A - 
(m 2 + 2 )/(u A lu 2 + 4) (Paczynski 1986; Griest 1991; Wambsganss 
2006), where u is the ratio of the AGN/lensing-galaxy separation 
(07 6) to the lensing galaxy’s Einstein radius, 
r 


Re = 


4GM\ ens Mens-lo- AGN 
C 2 Mens^AGN 


171 (Mens 12) 1/2 , 


where G is the gravitational constant, M is the lensing galaxy 
mass, M ens = 1300 Mpc is the angular size distance to the lens, 
Aagn - Da - 1400 Mpc is the angular size distance to the 
AGN, Mens-lo-AGN — MaGN “ (1 3 Z| cn s)/( 1 + ZaGN (M ens IS the 


distance between the lens located at redshift zi ens and the AGN 
located at redshift z A gn, and Mens 12 is M expressed in the units 
of 10 12 M o . The above amplification factor estimate involves a 
number of simplifications including (i) the simplified lens ge- 
ometry, (ii) use of the observed AGN-lens separation which is 
larger than the true one, and (Hi) the Paczynski (1986) formula 
referring to the combined light of two images (we know from 
observations that the single observed image of GB 1310+487 is 
much brighter than its second undetected image). Taking into ac- 
count these caveats, we estimate that the AGN image is probably 
amplified by a factor of a few. 

Since no second image is visible in optical Keck and radio 
VLBA (this work) and VLA (JVAS survey; King et al. 1999) 
data, we assume that its contribution to the source lightcurve at 
other bands, including GeV, is also negligible. The absence of an 
observable second image may indicate that either the lens is not 
massive and the source is still outside its Re or the mass distri- 
bution in the lens is asymmetric and the lensed source is close to 
a fold or cusp caustic. If the lens is a singular isothermal sphere 
(SIS; e.g., Refsdal & Surdej 1994; Meylan et al. 2006) and the 
lensed AGN is just outside the Einstein radius defined for an SIS 
through the lensing galaxy’s velocity dispersion cr SIS , 


Re sis 


= 47T 


°~SIS ^AGN 

C 2 Mens 


then a single image is formed having, in principle, an arbitrarily 
large amplification factor A — 1/(1 - 1/m) (e.g., Wu 1994). Tak- 
ing sis ^ 076 we estimate cr SIS < 140 km s -1 and the mass 
inside R^sis of < 6 X 10 10 M o (Fort & Mellier 1994). The low 
lensing galaxy mass, necessary to put the AGN image outside Re 
(and form a single image), may be reconciled with its brightness 
if the galaxy is undergoing intensive star formation, as indicated 
by strong emission lines in its spectrum. 

In general, gravitational lenses producing a single magnified 
image of a distant source should be more common than lenses 
producing multiple images. However, most gravitational lens 
searches (like the JVAS-CLASS survey; Browne et al. 2003) are 
designed to identify only multiple-image lenses. The BL Lac- 
type object AO 0235 + 164 is an example of a blazar shining 
through an intervening galaxy and having a single image weakly 
amplified by macrolensing (Abraham et al. 1993). 

A possibility of microlensing by individual foreground stars 
in the z = 0.5 galaxy cannot be excluded. The timescale of such 
microlensing events may be estimated as a ratio of the Einstein 
radius for a single star to the proper motion of the lens and is on 
the order of tens of years. Therefore, microlensing is probably 
unrelated to the observed fast high-energy variability, but may 
provide a significant amplification that is nearly constant over 
the duration of our observations. A large constant amplification 
due to microlensing could also explain the absence of the second 
lensed image. 


4.8. Emission model constraints from the SED 

The high-energy SED hump dominates over the synchrotron 
hump during the first (brighter) flare by a factor of q > 10. This 
is commonly observed in FSRQ-type blazars. In the framework 
of the leptonic model, the large q suggests that most of the ob- 
served y - ray flux during the first flare should be attributed to 
the EC process rather than to the SSC scenario. An SSC model 
would require a large deviation from equipartition: the emitting- 
particle energy density would need to be q 2 times greater than the 
magnetic field energy density in the source frame (Sikora et al. 
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2009), which is not expected; EC models do not require devia- 
tion from equipartition to explain large values of q. 

The dramatic difference in spectral slopes in the radio and IR 
regions suggests the presence of a break in the electron energy 
distribution. At the post-flare state, GB 1310+487 shows a steep 
y-ray spectrum. If this spectrum is dominated by the SSC emis- 
sion (Sect. 4.9), it is possible to estimate the Lorentz factor, y b , 
of electrons emitting at the energy distribution break from the 
positions of the SSC (vssc ~ 10 22 Hz) and synchrotron (v syn ~ 
10 14 Hz) SED peaks (Sikora et al. 2009): y b ~ (vssc/vsyn) 1/2 ~ 
10 4 . This value is an order of magnitude larger than those found 
in detailed SED modeling of FSRQs, (e.g., Vercellone et al. 
2011; Hayashida et al. 2012; Dutka et al. 2013). Different SED 
models applied to the same source may lead to different esti- 
mates ofy b (Sokolovsky et al. 2010; Abdo et al. 2010a). 

This is also the case for GB 1310+487 (Table 2). Hadronic 
models predict that X-rays are produced by synchrotron radi- 
ation of the secondary ultra-relativistic population of electrons 
and positrons. To reproduce the hard X-ray spectra observed in 
GB 1310+487 (Table 2) and in FSRQs, an extremely efficient ac- 
celeration of relativistic protons within the inner parts of the out- 
flow is needed, and the jet kinetic power must be orders of mag- 
nitude larger than the Eddington luminosity (Sikora et al. 2009), 
making this scenario very unlikely. 

4.9. Interpretation of changes in the SED 

According to the leptonic interpretation of blazar SEDs outlined 
in Sect. 1, the near-lR flux should be dominated by the syn- 
chrotron radiation while the observed y-ray flux is a combina- 
tion of EC and SSC components. The lack of broad lines in 
the optical spectrum observed by us (Sect. 3.4) and Shaw et al. 
(2012) suggests that the BLR in GB 1310+487 is weak or ob- 
scured from view. This, however, does not exclude the EC sce- 
nario: photons from the accretion disk or dusty torus might serve 
as targets for inverse-Compton scattering. The optical-UV spec- 
trum is flatter than the near IR one; however, due to contam- 
ination of the host (or intervening) galaxy and a nearby star 
(Sect. 3.3), it is not possible to distinguish accretion-disk emis- 
sion that might be present in this wavelength range. It may also 
be that the accretion disk avoids detection because it emits at 
shorter wavelengths. This would be the case if the central black 
hole mass is smaller than the one typically found in blazars, 
since the accretion-disk temperature decreases with increasing 
M. (Shakura & Sunyaev 1973). Far-IR data available only dur- 
ing the high-lR state (Flare 2) are not sufficient to estimate the 
possible contribution from a dusty torus. Thus, accretion-disk 
and dusty torus luminosities remain as free parameters in this 
discussion. 

The EC component peaks at higher energies than the SSC 
component for the following reasons. First, the typical energy 
of seed photons for the EC process (IR, optical, or UV corre- 
sponding to dusty toms, BLR, or accretion disk as the dominat- 
ing source of external radiation) should be higher than the typical 
energy of synchrotron photons. The synchrotron emission peak 
is located in the far-IR, as suggested by the observed steep near- 
IR spectrum. Second, due to relativistic aberration, most external 
photons illuminate the synchrotron-emitting plasma blob head 
on. Therefore, the external photons are additionally blueshifted 
in the reference frame of the plasma blob. A change in relative 
strength of the EC and SSC humps may explain the observed 
y-ray spectrum evolution, with the harder spectrum correspond- 
ing to greater contribution of the EC component to the total GeV 


flux. Large-amplitude GeV variability makes the corresponding 
spectral changes apparent even with the limited photon statistics 
of the observations. 

Another way to explain the observed changes in the y-ray 
spectrum is a varying contribution from multiple EC components 
(e.g., EC on accretion disk and dusty torus photons). This sce- 
nario, however, gives no predictions about the behavior of the 
synchrotron SED component, while the SSC+EC explanation is 
able to describe qualitatively the observed changes in the low- 
energy SED hump. 

The difference in SED during Flare 1 and Flare 2 may be 
qualitatively understood if the two flaring events are triggered 
by different physical mechanisms. The inverse-Compton hump 
brightening not associated with brightening of the synchrotron 
emission observed in Flare 1 may result from an increasing am- 
bient photon field that might be caused by an increasing accre- 
tion rate onto the central supermassive black hole (Paggi et al. 
201 1). 16 Flux increase in both synchrotron and inverse-Compton 
SED components (as observed in GB 1310+487 during Flare 2), 
combined with a peak energy increase of the two components, 
may result from additional electron acceleration. The inverse- 
Compton peak energy increase (with respect to pre- and post- 
flare states) is evident during Flare 2. The synchrotron peak en- 
ergy increase during Flare 2 is not excluded by the available data. 

This fits the pattern of y-ray spectrum changes discussed 
in Sect. 4.1, if the observed y-ray emission is a combination 
of EC emission peaking at higher energies and SSC emission 
peaking at lower energies. The first flare leads to the increased 
EC flux (relative to the SSC flux), which makes the overall y- 
ray spectrum harder. The second flare, probably caused by ad- 
ditional electron acceleration, is characterized by the enhanced 
synchrotron (as observed in the IR) and the corresponding SSC 
flux, possibly together with the EC component. The SSC com- 
ponent brightening makes the overall y-ray spectrum softer com- 
pared to Flare 1 . The post-flare high-energy hump might be a 
combination of SSC and EC components peaking at lower ener- 
gies than during Flare 2 due to lower acceleration of the under- 
lying electron population. Alternatively, since during the post- 
flare state the Compton dominance is less than an order of mag- 
nitude, the inverse-Compton hump might be dominated by the 
SSC component; SSC emission with q > 1 (as in Flare 2) might 
be produced (Zacharias & Schlickeiser 2012a,b) if the energy 
density of emitting particles is larger than that of the magnetic 
field, i.e., jet plasma is not in equipartition and Compton losses 
dominate the particle energy loss budget (e.g., Potter & Cotter 
2013). A constant supply of energy to the emitting particles is 
needed for this condition to be fulfilled for an extended period 
of time (Readhead 1994). Coordinated brightening of both syn- 
chrotron and Compton emission components during Flare 2 sug- 
gests that the SSC mechanism is the one responsible for the en- 
hanced Compton emission. 

It is possible that the lower-energy Compton emission is an- 
other EC component, not the SSC component. For example, the 
high-energy EC component responsible for Flare 1 might be as- 
sociated with accretion-disk photons, while the lower-energy EC 
component could be associated with lower-energy photons pro- 
duced in the dusty torus. Since the SSC emission has the same 

16 The increased accretion rate should result in accretion-disk brighten- 
ing. However, we may still not detect the accretion disk for the same 
reasons we do not detect it in the quiescent state: it is either too faint or 
emits at shorter wavelengths. Finally, we cannot exclude the possibility 
that the single available UV data point in the M2 band obtained during 
Flare 1 (Table 2, Fig. 10) has no contribution from the accretion-disk 
emission. 
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beaming pattern as synchrotron emission (e.g., Finke 2013), the 
SSC Compton dominance is independent of 6. The observed EC 
radiation intensity has a stronger dependence on 6 than does 
the synchrotron radiation (e.g., Georganopoulos et al. 2001). An 
increase in 6 would enhance both synchrotron and EC radi- 
ation, and the EC emission would be enhanced by a larger 
factor. This would explain the large value of q during Flare 2 
without the necessity of resorting to nonstationary SSC mod- 
els. The SSC peak in this scenario is assumed to be at even 
lower energies than the low-energy EC peak to account for the 
spectral change between Flare 2 and pre/post-flare states. This 
“SSC+low-energy-EC+high-energy-EC" scenario is more com- 
plex than the “SSC+EC" scenario, so we consider the former to 
be less likely. Additional detailed SED modeling is needed to 
test the above scenarios numerically. 

The intermediate interflare state is probably a combination of 
the decaying Flare 1 and rising Flare 2. This may indicate that ei- 
ther two independent emission zones are responsible for the two 
flares, or that a single emission region is propagating through a 
medium, with properties gradually changing from those which 
caused Flare 1 to those corresponding to Flare 2. 

The fact that Flare 1 has no counterpart in the 15 GFIz radio 
lightcurve (if a typical value of a few months for the y-ray/radio 
delay is assumed; see Sect. 4.5) suggests that the region responsi- 
ble for Flare 1 is located upstream of the VLBI core region. This 
region is probably heavily self-absorbed and does not contribute 
to the observed flux density at this frequency. Flare 2 occurred 
during the rise of the major radio flare; this may be an indication 
that the flaring region that dominates IR and high-energy emis- 
sion may be located close to the 15 GFIz radio core, as suggested 
for other blazars (e.g., Jorstad et al. 2010; Schinzel et al. 2012; 
Wehrle et al. 2012). 

5. Conclusions 

We summarize the major conclusions of this study as follows. 

1. Identification of the flaring y-ray source with the radio 
source GB 1310+487 is firmly established through positional 
coincidence and simultaneous multiwavelength observations 
of the flux variability. 

2. Significant changes in the y- ray photon index with flux were 
observed, showing the harder- when-brighter trend. It may re- 
sult from a changing relative contribution of EC and SSC 
emission to the total y- ray flux in the Fermi / LAT band. 

3. The bright near-IR flare does not correspond to the brightest 
y- ray flare. This may be an indication that different mecha- 
nisms are driving the two observed flares. 

4. From the absence of VLBA proper motion and the y - y 
opacity argument, we constrain the source Doppler factor: 
1.5 < 8 < 11. 

5. No clear association could be established between the y- ray 
variability and changes in radio flux and parsec-scale struc- 
ture. Simultaneous radio/y-ray observations over a longer 
time baseline are needed to test the possibility that some y- 
ray events are associated with radio flares. 

6. In the optical band, the object is a blend of a y-ray/radio-loud 
narrow-line AGN at z = 0.638 with an unrelated emission- 
line galaxy at z = 0.500. The AGN is not a member of com- 
mon types of y- ray flaring AGNs (blazars and NLSyls). 

7. The AGN radiation is probably amplified by a factor of a few 
because of gravitational lensing. 
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